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THE DYNAMICS OF WEAKLY REVERSIBLE POPULATION 
PROCESSES NEAR FACETS* 

DAVID F. ANDERSONt AND ANNE SHIU* 

Abstract. This paper concerns the dynamical behavior of weakly reversible, deterministically 
modeled population processes near the facets (codimension-one faces) of their invariant manifolds and 
proves that the facets of such systems are "repelling." It has been conjectured that any population 
process whose network graph is weakly reversible (has strongly connected components) is persistent. 
We prove this conjecture to be true for the subclass of weakly reversible systems for which only facets 
of the invariant manifold are associated with semilocking sets, or siphons. An important application 
of this work pertains to chemical reaction systems that are complex-balancing. For these systems it is 
known that within the interior of each invariant manifold there is a unique equilibrium. The Global 
Attractor Conjecture states that each of these equilibria is globally asymptotically stable relative to 
the interior of the invariant manifold in which it lies. Our results pertaining to weakly reversible 
systems imply that this conjecture holds for all complex- balancing systems whose boundary equilibria 
lie in the relative interior of the boundary facets. As a corollary, we show that the Global Attractor 
Conjecture holds for those systems for which the associated invariant manifolds are two-dimensional. 

Keywords: persistence, global stability, dynamical systems, population processes, chemical reaction 
systems, mass action kinetics, deficiency, complex-balancing, detailed-balancing, polyhedron. 

1. Introduction. Population processes are mathematical models that describe 
the time evolution of the abundances of interacting "species." To name a few exam- 
ples, population processes can be used to describe the dynamics of animal populations, 
the spread of infections, and the evolution of chemical systems. In these examples, 
the constituent species are the following: types of animals, infected and non-infected 
individuals, and chemical reactants and products, respectively. How best to model 
the dynamics of a population process depends upon the abundances of the constituent 
species. If the abundances are low, then the randomness of the interactions among 
the individual species is crucial to the system dynamics, so the process is most appro- 
priately modeled stochastically. On the other hand, if the abundances are sufficiently 
high so that the randomness is averaged out at the scale of concentrations, then the 
dynamics of the concentrations can be modeled deterministically. For precise state- 
ments regarding the relationship between the two models, see [171 118) . In the present 
paper we consider deterministic models. Also, we shall adopt the language associated 
with (bio)chemical reaction systems, which form a class of dynamical systems that 
arise in systems biology, and simply note that our results apply to any population 
process that satisfies our basic assumptions. 

The present work builds upon the body of work (usually called "chemical reac- 
tion network theory" ) that focuses on the qualitative properties of chemical reaction 
systems and, in particular, those properties that are independent of the values of the 
system parameters. See for example [3 [9l [131 US] . Examples of reaction systems from 
biology include pharmacological models of drug interaction [T^] , T-cell signal trans- 
duction models [SJ Uni [25] , and enzymatic mechanisms [22 • This line of research is 
important because there are many biochemical reaction systems that may warrant 
study at one time or another, and these systems are typically complex and highly 
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nonlinear. Further, the exact values of the system parameters are often unknown, 
and, worse still, these parameter values may vary from cell to cell. However, in a 
way that will be made precise in Section [21 the network structure of a given sys- 
tem induces differential equations that govern its dynamics, and it is this association 
between network structure and dynamics that can be utilized without the need for 
detailed knowledge of parameter values. 

To introduce our main results, we recall three terms from the literature that 
will be defined more precisely later. First, a directed graph is said to be weakly 
reversible if each of its connected components is strongly connected. The directed 
graphs we consider in this paper are chemical reaction diagrams in which the arrows 
denote possible reactions and the nodes are linear combinations of the species which 
represent the sources and products of the reactions. Second, for the systems in this 
paper a given trajectory is confined to an invariant polyhedron, which we shall denote 
by P. Such a polyhedron is called a positive stoichiometric compatibility class in the 
chemical reaction network theory literature and the faces of its boundary are contained 
in the boundary of the positive orthant. Third, semilocking sets, or siphons in the 
Petri net literature [21 HI [H] , are subsets of the set of species that characterize which 
faces of the boundary of P allow for the existence of equilibria and w-limit points. 

The main result of this paper, Theorem l3.21 concerns the dynamics of weakly re- 
versible chemical reaction systems near facets of P; a facet is a codimension-one face 
of P. Informally, Theorem 13.21 states that weak reversibility of the reaction diagram 
guarantees the following: for each point z found within the interior of a facet of P, 
there exists an open ( relative to P) neighborhood of z within which trajectories are 
forced away from the facet. Thus, Theorem 13.21 shows that weak reversibility guaran- 
tees that all facets are "repelling." We will prove this theorem by demonstrating that 
for each facet there must exist a reaction that pushes the trajectory away from that 
facet and that the corresponding reaction rate dominates all others. 

The main qualitative results of this paper concern the long term behavior of 
systems, and as such we are interested in the set of cj-limit points (accumulation 
points of trajectories). A bounded trajectory of a dynamical system for which R^q is 
forward invariant is said to be persistent if no tj-limit point lies on the boundary of 
the positive orthant. Thus, persistence corresponds to a non-extinction requirement. 
It has been conjectured that weak reversibility of a chemical reaction network implies 
that trajectories are persistent (for example, see [7]). Theorem 13.21 allows us to prove 
our main qualitative result. Theorem 13.41 which shows that this conjecture is true 
for the subclass of weakly reversible systems for which only facets of the invariant 
manifold are associated with semilocking sets. We also point out in Corollarv l3.5l that 
a slight variant of our proof of Theorem 13.21 shows that semilocking sets associated 
with facets are "dynamically non-emptiable" in the terminology of D. Angeli et al. 
[5], thereby providing a large class of dynamically non-emptiable semilocking sets. 

An important application of our main results pertains to chemical reaction sys- 
tems that are detailed-balancing or, more generally, complex-balancing [16]; these terms 
will be defined in Section^ For such systems, it is known that there is a unique equi- 
librium within the interior of each positive stoichiometric compatibility class P. This 
equilibrium is called the Birch point in [6] due to the connection to Birch's Theorem 
in Algebraic Statistics '201 Section 2.1]. Moreover, a strict Lyapunov function exists 
for this point, so local asymptotic stability relative to P is guaranteed [H [Ej. An 
open question is whether all trajectories with initial condition in the interior of P 
converge to the unique Birch point of P. The assertion that the answer is 'yes' is the 
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content of the Global Attractor Conjecture [SI [TB], 

The Global Attractor Conjecture is a special case of the conjecture discussed 
earlier that pertains to weakly reversible systems; in other words, one must show 
that all complex-balancing systems are persistent [7]. It is known that the set of 
Lj-limit points of such systems is contained within the set of equilibria [SI [53], so 
the conjecture is equivalent to the statement that any equilibrium on the boundary 
of P is not an w-limit point of an interior trajectory. Recent work has shown that 
certain boundary equilibria are not w-limit points of interior trajectories. For example, 
vertices of a positive stoichiometric compatibility class P are not w-limit points of 
interior trajectories even if they are equilibria [21[B]. In addition, the Global Attractor 
Conjecture recently has been shown to hold in the case that the system is detailed- 
balancing, P is two-dimensional, and the system is conservative (meaning that P is 
bounded) [6 . It is known that the underlying network of any detailed- or complex- 
balancing system necessarily is weakly reversible. Therefore, all of the results of 
this paper apply in this setting and when combined with previous results [5', W , give 
our main contribution to the Global Attractor Conjecture, Theorem 14.61 the Global 
Attractor Conjecture holds for systems for which the boundary equilibria are confined 
to facet-interior points and vertices of P. As a direct corollary to Theorem 14.61 we 
can conclude that the Global Attractor Conjecture holds for all systems for which P 
is two-dimensional, in other words, a polygon, thereby extending the result in [5]. 

We now describe the layout of the paper. Section [2| develops the mathematical 
model used throughout this paper. In so doing we also present concepts from poly- 
hedral geometry (Section [23]) that will be useful to us and formally define the notion 
of persistence (Section 12. 4p . In addition, the concept of a semilocking set is recalled. 
Our main results are then stated and proven in Section [3| Applications of this work to 
the Global Attractor Conjecture is the topic of Section [H Finally, Section [5| provides 
examples that illustrate our results within the context of related results. 

2. Mathematical formulation. In Sections 1 2. II and [2T2| we develop the mathe- 
matical model used in this paper and provide a brief introduction to chemical reaction 
network theory. In Section [2.3[ we present useful concepts from polyhedral geometry. 
In Section [2^ we recall the notions of persistence and semilocking sets. Throughout 
the following sections, we adopt the notation [n] := {1, 2, . . . , n}, for positive integers 
n G Z>o. 

2.1. Chemical reaction netviforks and basic terminology. An example of 
a chemical reaction is denoted by the following: 

2Xi + — > X2 ■ 

The Xi are called chemical species and 2Xi+Xt, and X2 are called chemical complexes. 
Assigning the source (or reactant) complex 2X1+X3 to the vector y = (2, 0, 3) and the 
product complex X2 to the vector y' = (0, 1, 0), we can write the reaction as y —>■ y' . 
In general we will denote by A'' the number of species Xi , and we consider a set of R 
reactions, each denoted by 

Vk y'k , 

for k G [R], and vectors yk,yk G ^>o; "^i^h yk y'k- Note that if = or y^ ~ 0, 
then this reaction represents an input or output to the system. Note that any complex 
may appear as both a source complex and a product complex in the system. For ease 
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of notation, when there is no need for enumeration we typically will drop the subscript 
k from the notation for the complexes and reactions. 

Definition 2.1. Let S = {Xi}, C = {y}, and TZ = {y ^ y'} denote sets 
of species, complexes, and reactions, respectively. The triple {S,C,Ti} is called a 
chemical reaction network. 

To each reaction network, {S,C,TZ}, we assign a unique directed graph (called a 
reaction diagram) constructed in the following manner. The nodes of the graph are 
the complexes, C. A directed edge (y, y') exists if and only if y ^ y' is a reaction in 
TZ. Each connected component of the resulting graph is termed a linkage class of the 
reaction diagram. 

Definition 2.2. The chemical reaction network is said to be weakly reversible 
if each linkage class of the corresponding reaction diagram is strongly connected. A 
network is said to be reversible if y' y G TZ whenever y ^ y' E TZ. Later we will 
say that a chemical reaction system is weakly reversible if its underlying network is. 

Let x{t) S denote the concentration vector of the species at time t with 
initial condition a;(0) = We will show in Section [^T^ that the vector x{t) — 
remains within the span of the reaction vectors {y'f. — yk}, i-e. in the linear space 
S = span{j/^ — y/c}feG[-R]> for time. We therefore make the following definition. 

Definition 2.3. The stoichiometric subspace of a network is the linear space 
S = span{yl - yk}ke[R]- 

It is known that under mild conditions on the rate functions of a system (see 
Section l^?^ . a trajectory x{t) with strictly positive initial condition x'^ G M^g remains 
in the strictly positive orthant M^q for all time (see Lemma 2.1 of 23 ). Thus, the 
trajectory remains in the open set {x'^ + 3)0 M^q, where x'^ + S := {z € M.^ \ z = 
x^ + V, for some v G S*}, for all time. In other words, this set is forward-invariant 
with respect to the dynamics. We shall refer to the closure of (x^ + S) DW^q, namely 

P (x" + 5)nM^o , (2.1) 

as a positive stoichiometric compatibility class. We note that this notation is slightly 
nonstandard, as in previous literature it was the interior of P that was termed the 
positive stoichiometric compatibility class. In the next section, we will see that P is 
a polyhedron. 

Remark 1. Ln spite of the notation, P clearly depends upon a choice of x^ £ ^>o- 
Throughout the paper, a reference to P assumes the existence of a positive initial 
condition x''^ G R^q f^''' which P is defined by (|2.ip . 

It will be convenient to view the set of species S as interchangeable with the set 
[N], where N denotes the number of species. Therefore, a subset of the species, W C 
iS, is also a subset of [N], and we will refer to the W -coordinates of a concentration 
vector X G M^, meaning the concentrations Xi for species i in W. Further, we will 
write i £ W or i £ [N] to represent Xi £ W or Xi £ S, respectively. Similarly, we 
sometimes will consider subsets of the set of reactions TZ as subsets of the set [R] . 

Definition 2.4. The zero-coordinates of a vector w £ are the indices i for 
which Wi — 0. The support of w is the set of indices for which Wi ^ 0. 

Based upon Definition 12.41 and the preceding remarks, both the set of zero- 
coordinates and the support of a vector w can, and will, be viewed as subsets of 
the species. 

2.2. The dynamics of a reaction system. A chemical reaction network gives 
rise to a dynamical system by way of a rate function for each reaction. In other words, 
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for each reaction y}~ we suppose the existence of a continuously differentiable 

function Rk{-) = Ry^^y'^{-) that satisfies the following assumption. 

Assumption 2.5. For k e [R], Rk{-) = Ryk^y'^{-) ■ K>o R satisfies: 

1. Ryi,^y' (■) depends explicitly upon Xi only if yki ^ 0. 

d 

2. — — Ryf,—,y>^{x) > for those Xi for which yki ^ 0, and equality can hold only 
ifx e dR^„. 

3. Ry^^y'Jx) = Oifxi = for some i with ykij^O- 

4- If ^ ^ Uki < yii, then lim \ — 0, where all other Xj > are held fixed 

xi-^o Rk(x) 

in the limit. 

The final assumption simply states that if the Zth reaction demands strictly more 
molecules of species Xi as inputs than does the fcth reaction, then the rate of the Ith 
reaction decreases to zero faster than the fcth reaction, as Xi 0. The functions Rk 
are typically referred to as the kinetics of the system and the dynamics of the system 
are given by the following coupled set of nonlinear ordinary differential equations: 



x{t)= M^my'k-yk) . (2.2) 

feG[-R] 

Integrating (|2.2p yields 

x{t) =x°+y2 (I Rk{x{s))ds\ [y'^ - yk) . 

Therefore, x{t) — x^ remains in the stoichiometric subspace, S = span{?;^ — yk}ke[R]i 
for all time, confirming the assertion made in the previous section. 

The most common kinetics, and the choice we shall make throughout the remain- 
der of this paper, is that of mass action kinetics. A chemical reaction system is said 
to have mass action kinetics if all functions Rk take the following multiplicative form: 

Rk{x) = Kkxf'xf^ ■ --xl'" =: tikx^' , (2.3) 

for some positive reaction rate constants Kk, where we have adopted the convention 
that 0° = 1 and the final equality is a definition. It is easily verified that each 
Rk defined via (j2.3|) satisfies Assumption 12.51 Combining (|2.2p and (|2.3p gives the 
following system of differential equations: 

x{t) = J2 ^'^^(t)"" iv'k - Vk) =■■ f{<t)) , (2.4) 

fcG[-R] 

where the last equality is a definition. This dynamical system is the main object of 
study in this paper. 

A concentration vector x G K>q is an equilibrium of the mass action system 
(|2.4p if f{x) = 0. Given that trajectories remain in their positive stoichiometric 
compatibility classes P for all positive time, we see that it is appropriate to ask about 
the existence and stability of equilibria of system (|2.4p within and relative to a positive 
stoichiometric compatibility class P. We will take this viewpoint in Section H) 

Remark 2. We note that every result in this paper holds for any chemical reaction 
systems with kinetics that satisfy Assumption \2.5[ Nonetheless, we choose to perform 
our analysis in the mass action case for clarity of exposition. 
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2.3. Connection to polyhedral geometry. We now recall terminology from 
polyhedral geometry that will be useful; we refer the reader to the text of G. Ziegler 
for further details I^. 

Definition 2.6. The half-space in M™ defined by a vector v e M™ and a constant 
c G M is the set 

Hy^c {xeM™ I {v,x)>c} . (2.5) 

A (convex) polyhedron in M™ is an intersection of finitely many half-spaces. 

For example, the non-negative orthant M>q is a polyhedron, as it can be written 
as the intersection of the N half-spaces He^^o, where the e^'s are the canonical unit 
vectors of R^. We now give three elementary facts about polyhedra from which we 
will deduce the fact that positive stoichiometric compatibility classes P are polyhedra. 
First, any linear space of R*" is a polyhedron. Second, any translation x -I- Q of a 
polyhedron Q by a vector x S K.™ is again a polyhedron. Third, the intersection of 
two polyhedra is a polyhedron. Therefore, as a translate -I- S) and the orthant 
R>Q are both polyhedra, it follows that the positive stoichiometric compatibility class 
P defined by (|2.ip is indeed a polyhedron. 

We continue with further definitions, which will allow us later to discuss boundary 
equilibria (those equilibria of (|2.4p on the boundary of P). 

Definition 2.7. Let Q be a polyhedron in R™. The interior ofQ, denoted by 
int((5), is the largest relatively open subset of Q. The dimension of Q, denoted by 
dim((5), is the dimension of the span of the translate of Q that contains the origin. 

For example, the dimension of P equals the dimension of the stoichiometric sub- 
space S: dim(P) = dim(5). We now define the faces of a polyhedron. 

Definition 2.8. Let Q be a polyhedron in R'". For a vector v S R™, the face 
of Q that it defines is the (possibly empty) set of points of Q that minimize the linear 
functional (w, ■) : R™ ^ R. 

If the minimum in Definition 12.81 (denoted Cmin) is attained, then we can write 
the face as F = Q n -ffw,c„,i„ H iJ_t,^c„i„- Therefore any face is itself a polyhedron, so 
we may speak of its dimension or its interior. 

Definition 2.9. Let Q be a polyhedron in R'". A facet of Q is a face whose 
dimension is one less than that of Q. A vertex is a nonempty zero- dimensional face 
(thus, it is a point). 

We make some remarks. First, note that what we call the "interior" is sometimes 
defined as the "relative interior" [S^. Second, vertices are called "extreme points" in 
[2] . Third, the interior of a vertex is seen to be the vertex itself. Fourth, the boundary 
of Q is the disjoint union of the interiors of the proper faces of Q. 

We now return to the positive stoichiometric classes P For a subset of the 

set of species W C S, let Zw C R^ denote its zero set: 

Zw = {xeR'^ : Xi = OifieW} . 

It can be seen that for any face F of a positive stoichiometric class P, there exists 
some possibly non-unique subset W d S such that 

F = Fw ■.= PnZw ■ (2.6) 

In other words, each face of P is the set of points of P whose set of zero-coordinates 
contains a certain subset W d S. However, it is important to note that for some sub- 
sets W, the face is empty: F\y — 0, and therefore no nonempty face of P corresponds 
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with such a W. In this case we say that the set Zw is stoichiometrically unattainable. 
We see also that Fw = P ii and only if W is empty. For definiteness, if there exist 
subsets Wi ^ W2 C S for which = Fw^ , we denote the face by Fw2 ■ Under this 
convention, it can be seen that the interior of a face Fw is 

int(Fiv) = { a; e P I a;, = if and only ifieW }. (2.7) 

We remark that the set mt{Fw) was denoted by Lw H P in [2]. 

The following example illustrates the above concepts. We note that in the interest 
of clarity we denote species by A, B,C, . . . rather than Xi,X2, X3, ... in all examples. 

Example 2.10. Consider the chemical reaction system which arises from the 
following reaction diagram: 

2A ^ A + B , B ^ C , (2.8) 

K2 K4 

where we use the standard notation of labeling a reaction arrow by the corresponding 
reaction rate constant. The stoichiometric subspace S in is spanned by the two 
reaction vectors (—1,1,0) and (0,-1,1). A positive stoichiometric compatibility class 
is depicted in Fiaure [2A\ it is a two-dimensional simplex (convex hull of three affinely 
independent points, in other words, a triangle) given by 

P= { {Xa,Xb,Xc) eM.io \ Xa+Xb+Xc--T } , (2.9) 

for positive total concentration T > 0. The three facets (edges) of each positive stoi- 
chiometric compatibility class P are one- dimensional line segments: 

F{A} = { (0, Xb, Xc) e M|o \xb + Xc^T} , 

F{B} = { (a^a, 0, Xc) e R|o \Xa+Xc^T] , 
F{C} = { {Xa,Xb, 0) e M>o \ Xa + Xb ^T } , 

and the three vertices are the three points P{a,b} = {(OiOj^)}^ P{A,c} — {(Oj^^jO)}; 
and F^B c} — Finally, the set Zs^a b,c} — {(05 0,0)} is stoichiometrically 

unattainable. We will revisit this reaction network in Examvle \5.1[ 
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2.4. Persistence and semilocking sets. Let x{t) be a solution to (|2.4p with 
strictly positive initial condition G K^g. The set of oj-limit points for this trajectory 
is the set of accumulation points: 

uj{x'^) := { X e M>o I x{tn) X, for some sequence t„ ^ oo }. (2-10) 

Definition 2.11. A bounded trajectory with initial condition x'^ is said to be 
persistent if uj{x^) H 9K>o = 0. A dynamical system with bounded trajectories is 
persistent if each trajectory with strictly positive initial condition is persistent. 

In order to show that a system is persistent, we must understand which points 
on the boundary of a positive stoichiometric class are capable of being cj-limit points. 
To this end, we recall the following definition from the literature. 

Definition 2.12. A nonempty subset W of the set of species is called a semilock- 
ing set if for each reaction in which there is an element ofW in the product complex, 
there is an element of W in the reactant complex. 

Remark 3. The notion of a semilocking set is the same as a siphon in the Petri 
net literature. See f3 \EIi- 

The intuition behind semilocking sets lies in the following proposition, which is 
the content of Proposition 5.5 of D. Angeli et al. [4]. Related results that concern 
the "reachability" of species include Theorems 1 and 2 in the textbook of A. Vol 'pert 
and S. Khudiaev [H Section 12.2.3]. 

Proposition 2.13 ([!]). Let W a S be non-empty. Then W is a semilocking 
set if and only if the face Fw is forward invariant for the dynamics ([23]). 

The above result holds because semilocking sets are characterized by the following 
property: if no species of W are present at time zero, then no species of W can be 
produced at any time in the future. In other words, these species are "locked" at 
zero for all time. If in addition the reaction network is weakly reversible, then it 
is straightforward to conclude the following: if a linkage class has a complex whose 
support contains an element of W , then the rates of all reactions within that linkage 
class will be zero for all positive time: = for any such reaction yk y'k- In 

other words, certain linkage classes are "shut off." 

In light of the characterization of the interior of a face Fw given in (|2.7p . the 
following theorem is proven in [2l|4]; it states that the semilocking sets are the possible 
sets of zero-coordinates of boundary cj-limit points. 

Theorem 2.14 (OH]). Let W C S be a nonempty subset of the set of species. 
Let x^ G be a strictly positive initial condition for the system (|2.4p . and let 

P ~ {x^ + S) n R'^q denote the corresponding positive stoichiometric compatibility 
class. If there exists an uj{x^) -limit point, z G uj^x'^), and a subset of the species, 
W , such that z is contained within the interior of the face Fw of P , then W is a 
semilocking set. 

Theorem 12 ■ 141 will be used in conjunction with results in the next section to prove 
the persistence of the following class of weakly reversible systems: those for which 
each semilocking set W satisfies dim(Fvv) = dim(P) — 1 (and so Fw is a facet of P) 
or Fw = (and so Zw is stoichiometrically unattainable); see Theorem 13.41 

3. Main results. In order to state Theorem l3.2i we need the following definition. 

Definition 3.1. Let Q c P be an open set relative to P, for which ^ ^ QC\dP CL 
Fw, for some face Fw of P. Then the face Fw is repelling in the neighborhood 
Q n int(P) with respect to the dynamics ()2.4I) if 

Y^xJ^ix) > (3.1) 
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for all X Q n int(P), where the fi are the functions given in Jg.^p . 

Remark 4. Note that Fyy is repelling in the neighborhood (5nint(P) with respect 
to if and only if ^ dist {x{t), Fw) > whenever x{t) G Q n int(F). Thus, Fw 

is repelling in a neighborhood if any trajectory in the neighborhood can not get closer 
to the face Fw while remaining in the neighborhood. 

The following theorem is our main technical result. 

Theorem 3.2. Let {S,C,TZ} be a weakly reversible chemical reaction network 
with dynamics governed by mass action kinetics (12. 4p . Let W d S be such that Fw 
is a facet of P , and take z to be in the interior of Fw ■ Then there exists a 5 > Q for 
which the facet Fw is repelling in the neighborhood Bs{z) nint(P), where Bs{z) is the 
open ball of radius 5 centered at the point z. 

Proof. For the time being we will assume that there is only one linkage class in the 
reaction diagram. The proof of the more general case is similar and will be discussed 
at the end. Also we assume that W — {Xi, . . . , Xm} for some M < N. 

Letting s dim(S') = dim(P), the facet Fw has dimension s — 1, which, 
by definition, means that Zw H 5 is an (s — l)-dimensional subspace of S. Let 
TT : M.^ be the projection onto the first M coordinates, that is, it is given by 

7r(xi, a;2, . . • , xm ,xm+i, ■ ■ ■ ,xn) ■— (a;i, X2, . . . , xm, 0, . . . , 0). As a shorthand we will 
also write x\w for Tr{x). Because Zw H S — ker(7r|s) has dimension (s — 1), it follows 
that the image 7t{S) is one-dimensional. Therefore, we may let u G S* be such that 
v\w spans the projection 7r(S'). We also let {w2, . . . ,Ws} span the subspace ZwC^S so 
that {v,W2, ... ,Ws} is SL basis for the subspace S. We simply note for future reference 
that by construction we have 



for each reaction y —> y' G TZ. Finally, for x e M.^ we define x|vf<= similarly to x\w', 
that is, x\w<: is the projection of x onto the final N — M coordinates. 

We may assume that all coordinates of are non-zero, for otherwise the concen- 
trations of certain species Xj G W would remain unchanged under the action of each 
reaction (note that we necessarily have Wfci = for fc G {2, . . . , s} and i G {1, . . . , M}). 
In such a case, the concentrations of those species Xj would remain constant in time, 
so we could simply remove them from the system by incorporating their concentrations 
into the rate constants appropriately. 

We will show that v\w has coordinates all of one sign and will use this fact to 
guarantee the existence of a "minimal complex" (with respect to the elements of 
W). We will then show that this minimal complex corresponds with a dominating 
monomial that appears as a positive term in each of the first M components of (|2.4p . 

Suppose, in order to find a contradiction, that v\w has coordinates of both positive 
and negative sign; that is, assume that Vi < < Vj for some indices i, j < M . Let 
u :— VjCi — ViCj G M>o (where ei denotes the Ith canonical basis vector). It follows 
that u G because (i) {u,v) = by construction, and (ii) {u,Wk) = for all 
k G {2, . . . , s} because these vectors have non-overlapping support. Note also that 
(w, z) = because the support of m is a subset of W whereas the support of z is W. 
Let a;° G R>g n P (such a point always exists by Remark [1]) . As z and both lie in 
P, there exist constants G M for fc G [s], such that 



y'\w - y\w e span(u|w) , 



(3.2) 



s 



z = X + aiv + 




k=2 
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Combining all of the above we conclude that 



s 

= {u, z) ~ {u,x^) + ai{u,v) + ''^^ak{u,Wk) = {u,x^) > 0, 

fe=2 



where the final inequality holds because u is non-negative and nonzero and x'^ has 
strictly positive components. This is a contradiction, so we conclude that v\w does 
not have both positive and negative coordinates, and, without loss of generality, we 
now assume that all coordinates of are positive. 

We recall from (|3.2p that y'\w — y\w G span(w|H') for each reaction y ^ y' G TZ. 
For concreteness we let y'i^\w — yk\w = "ykv\w for some 7fc G R where k S [R]. 
Combining this with the fact that Vi > for each i G {1, . . . , M} shows that each 
reaction yields either (i) a net gain of all species of W, {ii) a net loss of all species of 
W , or (iii) no change in any species of W and, moreover, that there exists a y G C 
such that y\w di y\w for all y E where we say x < y lox x,y & if Xi < yi 
for each i. Note that it is the sign of 7^. that determines whether a given reaction 
accounts for an increase or decrease in the abundances of the elements of W . 

We now find a neighborhood of positive radius 5 around z, Bs{z), for which the 
facet Fy/ is repelling in the neighborhood Bs{z) n int(P). The first condition we 
impose on S is that it must be less than the distance between z and any proper face 
of P that is not Fyy, which can be done because z is in the interior of the facet. 
Also, this condition ensures that for any point x £ Bs{z) n int(P), the coordinates 
Xi, for i > M, are uniformly bounded both above and below. Therefore, there exist 
constants i^min arid i^max such that for all x G Bs{z) D int(P) and all complexes y, 
we have the inequalities 

max ■ 

(3.3) 



The monomial x\^ will dominate all other monomials for x G Bs{z) D int(P) 
for sufficiently small S, and this will force trajectories away from the facet. To make 
this idea precise, let R+ denote those reactions that result in a net gain of the species 
in W and i?_ those that result in a net loss. We will write yk ^ y'k & R+ and 
yk ^ y'k ^ R- to enumerate over those reactions. We now have that for i G [M] and 
X G Bs{z) n int(P) and for sufficiently small S > 0, 



Mx) = V, ^^nkxfwxlpr - V, hkln.xf^^'^xf^ 

Vk^v'k<iR+ Vk^v'k^R- 



Finally, by weak reversibility (and possibly after choosing a different y that still sat- 
isfies the minimality condition), there is a reaction y y' E TZ with iji < y^ for all 
i G {!,... ,M} (i.e. the 7^ associated with this reaction is strictly positive). This 
reaction, which belongs to has a monomial, x\^^ , that necessarily dominates all 
monomials associated with reactions in i?_ (which necessarily have source complexes 
that contain a higher number of each element of W than y does). Combining this fact 
with (|3.4p shows that fi{x) > for i G [M] and for x G Bs{z) nint(P), and therefore, 
the facet Fw is repelling in the neighborhood Bs{z) n int(P). 
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In the case of more than one Hnkage class, each hnkage class will have its own 
minimal complex that will dominate all other monomials associated with that linkage 
class. Thus the desired result follows. □ 

Remark 5. Note that weak reversibility was used in the previous proof only to 
guarantee the existence of the reaction y ^ y' , where y\w is minimal and jji < y[ 
for all i G {1, . . . , Af }, and was not needed to prove the existence of such a y. If the 
network were not weakly reversible, but such a reaction nevertheless existed, then the 
same proof goes through unchanged. 

Corollary 3.3. Let {S,C,TZ} be a weakly reversible chemical reaction network 
with dynamics governed by mass action kinetics (|2.4p such that all trajectories are 
bounded. Suppose there exists a subset W CZ S , a positive initial condition x^^ G Il^>o, 
and a point z G u(xP^ fl Fy/ such that Fy/ is a facet of P. Then ti;(a;°) H dF^y ^ 0. 

Proof. Suppose not. That is, suppose that llj{x''^) n Fw C mt{Fw) holds. Let 
3^ :— w(a;°) n Fw- We claim that 3^ is a compact set; indeed, the trajectory x{t) is 
bounded so y is as well, and y is the intersection of two closed sets, and therefore is 
itself closed. 

Combining the compactness of 3^ C mt{Fw) with Theorem 13. 2[ we obtain that 
there exists an open covering of y consisting of a finite number of balls Bs.{zi) of 
positive radius Si, each centered around an element Zi of y, such that (z) each St 
is sufficiently small so that Bsi{zi) D dFy/ — 0, and (ii), for Q :— [JiBs^{zi), the 
facet Fw is repelling in Q n int(P). Combining these facts with the existence of 
z G uj{x°) n mt{Fw) n Q shows the existence of a point w G u;{x^) Ci int(Fw) fl dQ. 
However, this is impossible because w G u!{x'^) D int(Fvi/) = y necessitates that 
w € y C mt{Q). Thus, the claim is shown. □ 

We may now present our main qualitative result. 

Theorem 3.4. Let {S,C,TZ} be a weakly reversible chemical reaction network 
with dynamics governed by mass action kinetics (12. 4p such that all trajectories are 
bounded. Suppose that for each semilocking set W , the corresponding face Fyy is 
either a facet (dim(Fw) — dim{P) — 1) or is empty. Then the system is persistent. 

Proof. This is an immediate consequence of Theorem 12.141 and Corollary 13. 31 □ 

3.1. Connection to dynamically non-emptiable semilocking sets. In [3] 

the notion of a "dynamically non-emptiable" semilocking set is introduced. If for a 
given semilocking set W we define two sets 

C{W) I ^ a G : w 

= akiy'k — Vk) satisfies w\w d: ^ and 
Te{W) := { ^ a G M>o : a, < ea,, ^i,j G [R] such that yi\w ^ yj\w } , 

then W is said to be dynamically non-emptiable if C{W) H J-^{W) = {0} for some 
e > 0. Here the notation z < z' means that all coordinates satisfy the inequality 
Zi < z[ and z ^ z' means that furthermore at least one inequality is strict. Intuitively, 
this condition guarantees that all the concentrations of species of a semilocking set 
can not simultaneously decrease while preserving the necessary monomial dominance. 
D. Angeli et al. proved that if every semilocking set is dynamically non-emptiable, 
and if another condition holds (see ^ for details), then the system is persistent. 

We next prove that equation p.4p and a slight variant of the surrounding argu- 
ment can be used to show that any semilocking set W associated with a facet of a 
weakly reversible system is dynamically non-emptiable. We therefore have provided 
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a large set of examples of dynamically non-emptiable semilocking sets. 

Corollary 3.5. Let {S,C,TZ} he a weakly reversible chemical reaction network 
with dynamics governed by mass action kinetics (|2.4|) . Then any semilocking set 
associated with a facet is dynamically non-emptiable. 

Proof. As in the proof of Theorem 13.21 we may assume that there is one linkage 
class. The case of more than one linkage class is similar. Let W — {Xi, . . . , Xm} be 
a semilocking set such that Fw is a facet. Let ^ a G ^^{W) for some e > which 
may be made smaller as needed. Let v G ]R>o, with > if i G [M], be as in the 
proof of Theorem 13. 2( that is, y'\w — y\w G span(?;|vi/) for each reaction y ^ y' G TZ 
and y'klw — yk\w = Jkv\w for some 7fc G M where fc G [R]. Then, for w as in the 
definition of C{W), and i?+ and i?_ defined similarly as in the proof of Theorem 13. 2|, 
we have 

w\w = "^akiyklw ~ yk\w) = v\w ^ afc7fe - v\w ^ ctk\lk\ , 

and so for i G {1, ... , M}, 

Wi ^ Vi ^ aklk - Vi ^ Q!fc|7fe| . (3.5) 
Vk^y'k^R+ yk-*y'k'£R- 

Just as in the proof of Theorem 13. 2( there exists a reaction, the £th say, yi y'g ^ R+ 
such that yi di y for all y G C, and yi ^ yk if yk ^ y'k & Thus, for e > small 
enough, and by the definition of Te{W), the absolute value of the entire negative term 
in (|3.5p is less than or equal to Wmina£7£, where Wmin := minimi : i G [M]}. Therefore 
we see that w\w h 0. If ai = 0, then the preceding argument shows that ak = 
for yk ^ y'k ^ R- and so > for all i G [M] since a ^ 0; hence a ^ C{W). For 
the case ai ^ 0, it is clear that Wi > for all i G [M]. Thus, we again obtain that 
a ^ C{W) and the result is shown. □ 

4. The Global Attractor Conjecture. In this section, we use the results of 
the previous section to resolve some special cases of the Global Attractor Conjecture 
of chemical reaction network theory. In particular, the main result of this section, 
Theorem 14.61 establishes that the conjecture holds if all boundary equilibria are con- 
fined to the facets and vertices of the positive stoichiometric compatibility classes. 
That is, the conjecture holds if the faces associated with semilocking sets are facets, 
vertices, or are empty. 

4.1. Complex-balancing systems. The Global Attractor Conjecture is con- 
cerned with the asymptotic stability of so-called "complex-balancing" equilibria. Re- 
call that a concentration vector x G IR>q is an equilibrium of (j2.4p if the differential 
equations vanish at x: f(x) = 0. For each complex r] € C we will write {k \ yk = r]} 
and {k \ y'/. = r]} for the subsets of reactions A: G 7?. for which 77 is the source and 
product complex, respectively. In the following definition, it is understood that when 
summing over the reactions {k \ y'j^ — r]}, yk is used to represent the source complex 
of the given reaction. 

Definition 4.1. An equilib rium X G IR>q of (j2.4p is said to be complex-balancing 
if the following equality holds for each complex rj d C: 

i^kixr'^ - • 

{k I yk=ri} {k I y'^=v} 
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That is, X is a complex-balancing equilibrium if, at concentration x, the sum of reaction 
rates for reactions for which r/ is the source is equal to the sum of reaction rates for 
reactions for which r/ is the product. A complex-balancing system is a mass action 
system that admits a strictly positive complex-balancing equilibrium. 

In [6] , complex-balancing systems are called "toric dynamical systems" in order to 
highlight their inherent algebraic structure. Complex-balancing systems are automati- 
cally weakly reversible [8] . There are two important special cases of complex-balancing 
systems: the detailed-balancing systems and the zero deficiency systems. 

Definition 4.2. An equilibrium x G M>q of a reversible system with dynamics 
given by mass action (|2.4p is said to be detailed-balancing if for any pair of reversible 
reactions ^ with forward reaction rate Hk and backward rate k'j., the following 
equality holds: 

Kk{x)y' = Kki^'/" ■ 

That is, X is a detailed-balancing equilibrium if the forward rate of each reaction equals 
the reverse rate at concentration x. A detailed-balancing system is a reversible system 
with dynamics given by mass action (|2.4p that admits a strictly positive detailed- 
balancing equilibrium. 

Properties of detailed-balancing systems are described by Feinberg in [lOj and by 
A. Vol'pert and S. Khudiaev in [24l Section 12.3.3]. It is clear that detailed-balancing 
implies complex-balancing. 

Definition 4.3. For a chemical reaction network {S,C,TZ}, let n denote the 
number of complexes, I the number of linkage classes, and s the dimension of the 
stoichiometric subspace, S. The deficiency of the reaction network is the integer 
n — I — s. 

The deficiency of a reaction network is non-negative because it can be interpreted 
as either the dimension of a certain linear subspace [9] or the codimension of a certain 
ideal [B] . Note that the deficiency depends only on the reaction network or the reac- 
tion diagram. It is known that any weakly reversible dynamical system (j2.4p whose 
deficiency is zero is complex-balancing, and that this fact is independent of the choice 
of rate constants Kk 9 . On the other hand, a reaction diagram with a deficiency that 
is positive may give rise to a system that is both complex- and detailed-balancing, 
complex- but not detailed-balancing, or neither, depending on the values of the rate 
constants Kk P [51 [TUl [T5] . 

4.2. Qualitative behavior of complex-balancing systems. Much is known 
about the limiting behavior of complex-balancing systems. In the interior of each 
positive stoichiometric compatibility class P for such a system, there exists a unique 
equilibrium x, with strictly positive components, and this equilibrium is complex- 
balancing. As previously stated, x is called the Birch point due to the connection 
to Birch's Theorem (see Theorem 1.10 of [^). Note that a system was defined to 
be complex-balancing if at least one such equilibrium exists; we now are asserting 
that so long as at least one P contains a complex-balancing equilibrium, then they 
all do. As for the stability of the equilibrium within the interior of the corresponding 
P, a strict Lyapunov function exists for each such point. Hence local asymptotic 
stability relative to P is guaranteed; see Theorem 6A of [TC] and the Deficiency Zero 
Theorem of [Hj. The Global Attractor Conjecture states that this equilibrium of P 
is globally asymptotically stable relative to the interior of P [B]. In the following 
statement, a global attractor for a set is a point v* ^ V such that any trajectory 
v{t) with initial condition v'^ V converges to v* , in other words, lim v(t) = v*. 

t — ^oo 
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Global Attractor Conjecture For any complex-balancing system (|2.4p and any 
strictly positive initial condition , the Birch point x E P :— {x^ + S)r\M.^Q is a global 
attractor of the interior of the positive stoichiometric compatibility class, int(P). 

This conjecture first appeared in a paper of F. Horn [14] . and was given the name 
"Global Attractor Conjecture" by G. Craciun et al. [Sj. It is stated to be the main 
open question in the area of chemical reaction network theory by L. Adleman et al. 
[T]. In fact, M. Feinberg stated the more general conjecture that all weakly reversible 
systems are persistent Section 6.1]. To this end, G. Gnacadja proved that the class 
of networks of "reversible binding reactions" are persistent; these systems include 
non-complex-balancing ones [11] . 

We now describe known partial results regarding the Global Attractor Conjec- 
ture. By an interior trajectory we shall mean a solution x{t) to (|2.4p that begins 
at a strictly positive initial condition x'^ S IR>o- known that trajectories of 

complex-balancing systems converge to the set of equilibria; see Corollary 2.6.4 of [5] 
or Theorem 1 of [23]. Hence, the conjecture is equivalent to the following statement: 
for a complex-balancing system, any boundary equilibrium is not an uj -limit point of 
an interior trajectory. It clearly follows that if a positive stoichiometric compatibility 
class P has no boundary equilibria, then the Global Attractor Conjecture holds for 
this P. Thus, sufficient conditions for the non-existence of boundary equilibria are 
conditions for which the Global Attractor Conjecture holds (see Theorem 2.9 of [2]); 
a result of this type is Theorem 6.1 of L. Adleman et al. [1]. Recall that by Theo- 
rem 12.141 we know that the only faces Fy/ of a positive stoichiometric compatibility 
class P that may contain w-limit points in their interiors are those for which W is 
a semilocking set. In particular, if the set Zw is stoichiometrically unattainable for 
all semilocking sets W , then P has no boundary equilibria, and hence, the Global 
Attractor Conjecture holds for this P; see the main theorem of D. Angeli et al. [4]. 
Biological models in which the non-existence of boundary equilibria implies global 
convergence include the ligand-receptor-antagonist-trap model of G. Gnacadja et al. 
[T2] . the enzymatic mechanism of D. Siegel and D. MacLean [52], and the T-cell signal 
transduction model of T. McKeithan [TH] (the mathematical analysis appears in the 
work of E. Sontag [23] and in the Ph.D. thesis of M. Chavez [3 Section 7.1]). We 
remark that this type of argument first appeared in [7l Section 6.1]. 

The remaining case of the Global Attractor Conjecture, in which equilibria exist 
on the boundary of P, is still open. However, some progress has been made. For 
example, it already is known that vertices of P can not be w-limit points even if 
they are equilibria; see Theorem 3.7 in the work of D. Anderson [2] or Proposition 20 
of the work of G. Craciun et al. [6]. For another class of systems for which the 
Global Attractor Conjecture holds despite the presence of boundary equilibria, see 
Proposition 7.2.1 of the work of M. Chavez [5]. The hypotheses of this result are 
that the set of boundary equilibria in P is discrete, that each boundary equilibrium 
is hyperbolic with respect to P, and that a third, more technical condition holds. 
In addition, the Global Attractor Conjecture holds in the case that the network is 
detailed-balancing, P is two-dimensional, and the network is conservative (meaning 
that P is bounded); see Theorem 23 of [6]. In the next section, Corollarv l4 . 71 will allow 
us to eliminate the hypotheses "detailed-balancing" and "conservative" from the two- 
dimensional result; such partial results concerning the Global Attractor Conjecture are 
the next topic of this paper. 

4.3. Applications to complex-balancing systems. Theorem l4.6l is our main 
contribution to the Global Attractor Conjecture and generalizes the known results 
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described in Section [4.21 We first present a definition and a lemma. 

Definition 4.4. Suppose that {S,C,TZ} is a weakly reversible chemical reac- 
tion network, endowed with mass action kinetics, and W G S is a semilocking set. 
Then, the W^-reduced network is the chemical reaction network {S \ W,C,TZ} such 
that C and TZ are those complexes and reactions that do not involve a species from 
W . Equivalently, this is the subnetwork obtained by removing all linkage classes that 
contain a complex with a species from W in its support. The M^-reduced system is the 
W -reduced network endowed with the same rate constants as in the original system. 

As noted in comments following Proposition 12.131 for a weakly reversible system 
and any semilocking set W , either each complex in a given linkage class contains an 
element of W or each complex in that linkage class does not contain any element of 
W . Therefore, W^-reduced systems are themselves weakly reversible. Furthermore, it 
is easy to check that any ly-reduced system of a complex-balancing system is itself 
complex-balancing. 

Lemma 4.5. Consider a complex-balancing system. A face F\y of a stoichiometric 
compatibility class P contains an equilibrium in its interior if and only if W is a 
semilocking set. 

Proof. It is clear that if a face Fw contains an equilibrium in its interior, then 
IV is a semilocking set (see the discussion following Proposition 12.131 or the proof of 
Theorem 2.5 of 0). 

If IV is a semilocking set, then the TV-reduced system is complex-balancing and 
one of its invariant polyhedra is naturally identified with Fw- Therefore, this reduced 
system admits a complex-balancing equilibrium, which we denote by y. For i G W, 
let Zi be the component of y associated with species i. For i S W, let Zi = 0. So 
constructed, the concentration vector z is an equilibrium within the interior of ^V- D 

Theorem 4.6. The Global Attractor Conjecture holds for any complex-balancing 
(and in particular, detailed-balancing or weakly reversible zero deficiency) chemical 
reaction system whose boundary equilibria are confined to facet-interior points or ver- 
tices of the positive stoichiometric compatibility classes. Equivalently, if a face Fw is 
a facet, vertex, or an empty face whenever W is a semilocking set, then the Global 
Attractor Conjecture holds. 

Proof. The equivalence of the two statements in the theorem is a consequence of 
Lemma |4?5] As noted in the previous section, persistence is a necessary and sufficient 
condition for the Global Attractor Conjecture to hold. Further, by the results in [2] 
or [6], vertices may not be w-limit points. The remainder of the proof is similar to 
that of Theorem 13.41 and is omitted. □ 

We note that Theorem 14.61 and previous results in [51 [S] show that for equilibria 
z of complex-balancing systems that reside within the interiors of facets or of vertices 
of P, these faces are repelling near z. More precisely, there exists a relatively open 
set Q of P such that Q contains z and the face is repelling in Q n int(P). 

The following corollary resolves the Global Attractor Conjecture for systems of 
dimension two; note that the one-dimensional case is straightforward to prove. 

Corollary 4.7 (GAG for two-dimensional P). The Global Attractor Conjec- 
ture holds for all complex-balancing (and in particular, detailed-balancing or weakly 
reversible zero deficiency) chemical reaction systems whose positive stoichiometric 
compatibility classes are two-dimensional. 

Proof. This follows immediately from Theorem l4.6l as each face of a two-dimensional 
polytope must be either a facet or a vertex. □ 

In the next section, we provide examples that illustrate our results, as well as a 
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three-dimensional example for which our results do not apply. 

5. Examples. As discussed in the previous section, the Global Attractor Con- 
jecture previously has been shown to hold if P has no boundary equilibria or if the 
boundary equilibria are restricted to vertices of P. Therefore, the examples in this 
section that pertain to complex-balancing systems feature non- vertex boundary equi- 
libria. Also, we have chosen examples for which the conditions of the theorems can 
be easily checked. 

Example 5.1. We revisit the network given by the following reactions: 
2A ^ A + B , B ^ C . 

As we saw in Examvle \2. 101 the positive stoichiometric compatibility classes are two- 
dimensional triangles: 

F= { {Xa,Xb,Xc) e Rio \ ^a+Xb+Xc^T } , (5.1) 

where T > 0. It is straightforward to check that P has a unique boundary equilibrium 
given by 

z = (o, -^T, -^T 

\ H3 + Hi K3 + K4 

and that this point lies in the interior of the facet ^{a}- (Note that this boundary 
equilibrium is the Birch point of the reversible deficiency zero subnetwork B ^ C.) 
Therefore, both Theorem \4-f}] and Corollary \4. 7| allow us to conclude that despite the 
presence of the boundary equilibrium z, the Birch point in the interior of P is globally 
asymptotically stable. 

We remark that the results in [2] do not apply to the previous example, although 
Theorem 23 of [6 and Theorem 4 of do. However, for the following example, no 
previously known results apply. 

Example 5.2. Consider the reaction network depicted below 

2A ^ A + B 



A + C 

The positive stoichiometric compatibility classes are the same triangles Ii5.1\) as in the 
previous example. For each P, the set of boundary equilibria is the entire face F^j^y 
(one of the three edges of P), which includes the two vertices -F{a,b} o,n-d ^{a,c}- 
Hence the results 0/ rfo not apply. Note that this is a weakly reversible zero 
deficiency network, but it is not detailed-balancing; so the results of do not apply. 
Finally, the second condition of Theorem 4 of f^' is not satisfied by this example 
(there are nested "deadlocks" f^) and so that result also does not apply. However, 
both Theorem \4-.b\ and Corollary \4-T\ imply that the Global Attractor Conjecture holds 
for all choices of rate constants and for all P defined by this network, despite the 
presence of boundary equilibria. 

In the next example, the positive stoichiometric compatibility classes are three- 
dimensional. 

Example 5.3. The following zero deficiency network is obtained from Exam- 
ple 15. Jl by adding a reversible reaction: 

2A U A + B , B U C U D . 
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The positive stoichiometric compatibility classes are three-dimensional simplices (tetra- 
hedra): 

P = { {Xa,Xh, Xc, Xd) e Ryo \ Xa + Xb + Xc + Xd ^ T } , 

for positive total concentration T > 0. The unique boundary equilibrium in P is the 
Birch point of the zero deficiency subnetwork B ^ C ^ D, and it lies in the interior 
of the facet ^{a}- In- other words, the point is z = {0,Xb,Xc,Xd) where {xi,,Xc,Xd) is 
the Birch point for the system defined by the subnetwork 

B C ^ D . 

Thus by Theorem \4-(>\ the Global Attractor Conjecture holds for all positive stoichio- 
metric compatibility classes P and all choices of rate constants defined by this net- 
work. 

As in the previous example, the positive stoichiometric compatibihty classes of 
our next example are three-dimensional. However neither previously known results 
[21 [5l [6] nor our current results can resolve the question of global asymptotic stability. 

Example 5.4. The following zero deficiency network consists of three reversible 
reactions: 

A B ^ A + B A + C . 

As there are no conservation relations, the unique positive stoichiometric compatibility 
class is the entire non-negative orthant: 



The set of boundary equilibria is the one- dimensional face (ray) -Fj^.s}; which includes 
the origin F{a,_b,c}- Therefore non-vertex, non-facet boundary equilibria exist, so the 
results in this paper do not apply. 

We end with an example in which the results of Section [3] apply but those of 
Section [4] do not. 

Example 5.5. The following reversible network is obtained from Example \5.1\ by 
adding another reversible reaction: 

A + C ^ 2A ^ A + B , B T± C . 

The positive stoichiometric compatibility classes are again the two-dimensional trian- 
gles given by (|5.ip . One can easily check that the network has a deficiency of one, so 
there exist rate constants for which the system is not complex-balancing (for example, 
Ki = K3, K5 = Kg, and K2 K4). Thus the results of Section^ do not apply. It is also 
easy to verify that {A} and {A,B,C} are the only semilocking sets and that ^{yi} 
a facet a?jrf F{a,s,c} empty. Therefore, Theorem \3.4\ applies and we conclude that, 
independent of the choice of rate constants, the system is persistent. 
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